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Abstract. In this paper we consider applications of methods from wavelet analysis 
Z_) • to nonlinear dynamical problems related to accelerator physics. In our approach we 

Cw I take into account underlying algebraical, geometrical and topological structures of 

C/2 ' corresponding problems. 

^'- I INTRODUCTION 

This paper is the sequel of our first paper in this volume [1], in which we con- 
sidered the applications of a number of analytical methods from nonlinear (local) 
^ • Fourier analysis, or wavelet analysis, to nonlinear accelerator physics problems. 

^ . This paper is the continuation of results from [2]-[7], which is based on our approach 

^ I to investigation of nonlinear problems both general and with additional structures 

O I (Hamiltonian, symplectic or quasi comp lex) , chaotic, quasiclassical, quantum. 

Wavelet analysis is a relatively novel set of mathematical methods, which gives 

us a possibility to work with well-localized bases in functional spaces and with the 

O I general type of operators (differential, integral, pseudodifferential) in such bases. 

c/3 ' In contrast with paper [1] in this paper we try to take into account before using 

power analytical approaches underlying algebraical, geometrical, topological struc- 

p |! tures related to kinematical, dynamical and hidden symmetry of physical problems. 

In this paper we give a review of a number of the corresponding problems and 

describe the key points of some possible methods by which we can find the full 

solutions of initial physical problem. We described a few concrete problems in [1, 

part II]. The most interesting case is the dynamics of spin-orbital motion [1, II D]. 

Related problems may be found in [8]. 

The content of this paper is not more than an attempt to extract the most 
complicated formal or mathematical or principal parts of the World of nonlinear 
accelerator physics, which is today beyond of mainstream in our opinion. 

In part II we consider dynamical consequences of covariance properties regarding 
to relativity (kinematical) groups and continuous wavelet transform as a method 



for the solution of dynamical problems. 

In part II A we introduce the semidirect product structure, which allows us to 
consider from general point of view all relativity groups such as Euclidean, Galilei, 
Poincare. 

Then in part II B we consider the Lie-Poisson equations and obtain the manifes- 
tation of semiproduct structure of (kinematic) symmetry group on dynamical level. 
So, correct description of dynamics is a consequence of correct understanding of 
real symmetry of the concrete problem. 

In part II C we consider the technique for simplifications of dynamics related 
to semiproduct structure by using reduction to corresponding orbit structure. As 
result we have simplified Lie-Poisson equations. 

In part II D we consider the Lagrangian theory related to semiproduct struc- 
ture and explicit form of variation principle and corresponding (semidirect) Euler- 
Poicare equations. 

In part II E we introduce a continuous wavelet transform and corresponding 
analytical technique which allow to consider covariant wavelet analysis. 

In part II F we consider in the particular case of afiine Galilei group with the 
semiproduct structure also the corresponding orbit technique for constructing dif- 
ferent types of invariant wavelet bases. 

In part III we consider instead of kinematical symmetry the dynamical symmetry. 

In part III A according to the orbit method and by using construction from the 
geometric quantization theory we construct the symplectic and Poisson structures 
associated with generalized wavelets by using metaplectic structure. We consider 
wavelet approach to the calculations of Melnikov functions in the theory of homo- 
clinic chaos in perturbed Hamiltonian systems in part III B and for calculation of 
Arnold- Weinstein curves (closed loops) in Floer variational approach in part III C. 

In parts III D, III E we consider applications of very useful fast wavelet transform 
technique (part III F) to calculations in symplectic scale of spaces and to quasi- 
classical evolution dynamics. This method gives maximally sparse representation 
of (differential) operator that allows us to take into account contribution from each 
level of resolution. 

In part IV A we consider symplectic and Lagrangian structures for the case of 
discretization of flows by corresponding maps and in part IV B construction of 
corresponding solutions by applications of generalized wavelet approach which is 
based on generalization of multiresolution analysis for the case of maps. 

II SEMIDIRECT PRODUCT, DYNAMICS, WAVELET 

REPRESENTATION 

A Semidirect Product 

Relativity groups such as Euclidean, Galilei or Poincare groups are the particular 
cases of semidirect product construction, which is very useful and simple general 



construction in the group theory [9]. We may consider as a basic example the 
Euchdean group SE{3) = 5*0(3) cxi R^, the semi direct product of rotations and 
translations. In general case we have S = G txi V, where group G (Lie group or 
automorphisms group) acts on a vector space V and on its dual V*. Let V he a 
vector space and G is the Lie group, which acts on the left by linear maps on V (G 
also acts on the left on its dual space V*). The semidirect product S = G (xV is 
the Cartesian product S = G x V with group multiplication 

i9i,vi)i92,V2) = i9i92,vi+ giV2), (1) 

where the action oi g & G on v & V is denoted as gv. Of course, we can consider 
the corresponding definitions both in case of the right actions and in case, when G 
is a group of automorphisms of the vector space V. As we shall explain below both 
cases. Lie groups and automorphisms groups, are important for us. 

So, the Lie algebra of S is the semidirect product Lie algebra, s = Q ix V with 
brackets 

m,V^),{^2,V2)] = ([ei,6],6^2-6^^l), (2) 

where the induced action of Q by concatenation is denoted as .^if2. Also we need 
expressions for adjoint and coadjoint actions for semidirect products. Let {g,v) G 
S = G xV, {^,u) e s = g xV, {fi,a) e s* = g* xV*, g^ = Adg^, gfi = Ad^-ifi, 
ga denotes the induced left action oi g on a (the left action of G on V induces a 
left action on V* — the inverse of the transpose of the action on V) , p^ : Q -^ V is 
a linear map given by Pv{0 — ^'^t Pl '■ ^* ^ G* is its dual. Then these actions are 
given by simple concatenation: 

{9,v){i,u) = {gi,gu-{gCjv), (3) 

{g,v){fi,a) = {gfi + pl{ga),ga) 

Below we use the following notation: p^a = v o a E Q* for a E V* , which is a 
bilinear operation in v and a. So, we have the coadjoint action: 

{g,v){p,a) = {gp + vo{ga),ga). (4) 

Using concatenation notation for Lie algebra actions we have alternative definition 
oi V o a E Q* . For dX\ v E V , a E V* , r] E Q we have 

< 77a, V >= — < V o a,!] > (5) 



B The Lie-Poisson Equations and Semiproduct Structure 

Below we consider the manifestation of semiproduct structure of symmetry group 
on dynamical level. Let F, G be real valued functions on the dual space Q* , p E Q* . 
Functional derivative of F at /i is the unique element SF/Sp G Q: 



1 5F 

lim -[F(/i + e(5/i) - F(/i)] =< (5/i, — > (6) 

for all (5/i G ^*, <, > is pairing between Q* and ^. 
Define the (±) Lie-Poisson brackets by 

{F,G}±(/i) = ±</i,[^,^]> (7) 

The Lie-Poisson equations, determined by 

F = {F,H} (8) 

read intrinsically 

A = Tad*gH/d^,f^- (9) 

For the left representation of G on V ± Lie-Poisson bracket of two functions /, k : 
s* ^ R is given by 

if, k}±{fM, a) = ± < /i, [—,—]> ± < a, — — - ^^ >, (10) 

ofi ojj, o^oa ojioa 

where Sf /6n G Q, 6f/6a G V are the functional derivatives of f (6). The Hamilto- 
nian vector field of /i : s* G R has the expression 

AA / N / ,* ^h Sh . .... 

Thus, Hamiltonian equations on the dual of a semidirect product are [9]: 

(x = ^ad*sy^ls^^± — oa (12) 

5h 
a = =b-— a 

OjJ, 

So, we can see the explicit difference between Poisson brackets (7) and (10) and 
the equations of motion (9) and (12), which come from the semiproduct structure. 

C Reduction of Dynamics on Semiproduct 

There is technique for reducing dynamics that is associated with the geometry 
of semidirect product reduction theorem[9]. Let us have a Hamiltonian on T*G 
that is invariant under the isotropy Ga^ for Oq ^ ^*- The semidirect product 
reduction theorem states that reduction of T*G by Gao gives reduced spaces that 
are simplectically diffeomorphic to coadjoint orbits in the dual of the Lie algebra 



of the semidirect product {Q ixi V)*. If one reduces the semidirect group product 
S = G ixi V in two stages, first by V and then by G one recovers this semidirect 
product reduction theorem. Thus, let S = G (xV, choose a = {/j,, a) ^ Q* x V* and 
reduce T*S by the action of S* at a giving the coadjoint orbit O^ through a & S*. 
There is a symplectic diffeomorphism between O^ and the reduced space obtained 
be reducing T*G by the subgroup Ga (the isotropy of G for its action on V* at the 
point a E V*) at the point nlQa, where Qa is the Lie algebra of Ga- 
Then we have the following procedure. 

1. We start with a Hamiltonian Ha^ on T*G that depends parametrically on a 
variable a^ E V*. 

2. The Hamiltonian regarded as a map: T*GxV* ^ R is assumed to be invariant 
on T*G under the action of G on T*^ x V*. 

3. The condition 2 is equivalent to the invariance of the function H defined on 
T*S = T*G X V X V* extended to be constant in the variable V under the 
action of the semidirect product. 

4. By the semidirect product reduction theorem, the dynamics of Hag reduced 
by Gao, the isotropy group of ao, is simplectically equivalent to Lie-Poisson 
dynamics on s* = Q* x V* . 

5. This Lie-Poisson dynamics is given by equations (12) for the function h{fi, a) = 
H{ag,g~^a), where fj, = g~^ag. 

D Lagrangian Theory, the Euler-Poincare Equations, 
Variational Approach on Semiproduct 

Now we consider according to [9] Lagrangian side of a theory. This approach is 
based on variational principles with symmetry and is not dependent on Hamilto- 
nian formulation, although it is demonstrated in [9] that this purely Lagrangian 
formulation is equivalent to the Hamiltonian formulation on duals of semidirect 
product (the corresponding Legendre transformation is a diffeomorphism). 

We consider the case of the left representation and the left invariant Lagrangians 
{£ and L), which depend in additional on another parameter a E V* (dynamical 
parameter), where V is representation space for the Lie group G and L has an 
invariance property related to both arguments. It should be noted that the result- 
ing equations of motion, the Euler-Poincare equations, are not the Euler-Poincare 
equations for the semidirect product Lie algebra ^ ixi \/* or ^ ixi K. 

So, we have the following: 

1. There is a left presentation of Lie group G on the vector space V and G acts 
in the natural way on the left on TG x V* : h{vg, a) = {hvg, ha). 

2. The function L : TG x \/* G R is the left G-invariant. 



3. Let ao G V*, Lagrangian La^ : TG -^ R, La^X'^g) = L{vg,aQ). Lag is left 
invariant under the lift to TG of the left action of Gao on G, where Gag is the 
isotropy group of ao- 

4. Left G-invariance of L permits us to define 

e:gxV*^R (13) 

by 

'^Vg, g'^ao) = L{vg, ao). (14) 



This relation defines for any i : Q x V* ^ H the left G-invariant function 
L:TGxV* ^R. 

5. For a curve g(t) EG let be 

m := 9{t)-'m (15) 

and define the curve a{t) as the unique solution of the following linear differ- 
ential equation with time dependent coefficients 

d(t) = -e(t)a(t), (16) 

with initial condition a(0) = ag. The solution can be written as a{t) = 
g{t)'^ao. 

Then we have four equivalent descriptions of the corresponding dynamics: 

1. If Co is fixed then Hamilton's variational principle 

6 f''LaMt),mdt = (17) 

Jti 

holds for variations Sg{t) of g{t) vanishing at the endpoints. 

2. g(t) satisfies the Euler-Lagrange equations for La^ on G. 

3. The constrained variational principle 

6 rmt),a{t))dt = (18) 

Jtl 

holds on ^ X V"*, using variations of ^ and a of the form 6$, = rj + [^iV]^ 
5a = —rja, where ?7(t) G Q vanishes at the endpoints. 

4. The Euler-Poincare equations hold on ^ x V"* 

— — - = ad.— : + -— oa (19) 

dt5^ H^ Sa ^ ' 

So, we may apply our wavelet methods either on the level of variational formulation 
(17) or on the level of Euler-Poincare equations (19). 



E Continuous Wavelet Transform 

Now we need take into account the Hamiltonian or Lagrangian structures related 
with systems (12) or (19). Therefore, we need to consider generahzed wavelets, 
which allow us to consider the corresponding structures instead of compactly sup- 
ported wavelet representation from paper [1]. 

In wavelet analysis the following three concepts are used now: 1). a square in- 
tegrable representation [/ of a group G, 2). coherent states (CS) over G, 3). the 
wavelet transform associated to U. We consider now their unification [10], [11]. 

Let G be a locally compact group and Ua strongly continuous, irreducible, unitary 
representation of G on Hilbert space Ti. Let iJ be a closed subgroup of G, X = G/ H 
with (quasi) invariant measure v and a : X = GjH ^ G is a Borel section in a 
principal bundle G -^ GjH . Then we say that \J is square integrable modiH^ cr) if 
there exists a non-zero vector rj & 7i such that 







< / I < U(a(x))ri\<^ > Pdz/(x) =< $1^^$ > < cx), V^ G H (20) 

Jx 



Given such a vector rj E 7i called admissible for {U, a) we define the family of 
(covariant) coherent states or wavelets, indexed by points x E X, as the orbit of r] 
under G, though the representation U and the section a [10], [11] 

S. = V.{.) = Uia{x))r]\x E X (21) 

So, coherent states or wavelets are simply the elements of the orbit under U of a fixed 
vector 1] in representation space. We have the following fundamental properties: 

1. Overcompleteness: 

The set S^ is total in H : (S'o-)"'" = 

2. Resolution property: 

the square integrability condition (20) may be represented as a resolution 
relation: 

/ \ri^{x) >< ri^(^^)\du{x) = A^, (22) 

where A^r is a bounded, positive operator with a densely defined inverse. Define 
the linear map 

W,:n^ L\X,du), {W,^){x) =< Vai.)\<^ > (23) 

Then the range H^ of Wrj is complete with respect to the scalar product 
< $|\l/ >jj=< ^\W^A~'^W~^'^ > and W^j is unitary operator from H onto 7i^. 
Wrj is Continuous Wavelet Transform (CWT). 



3. Reproducing kernel 

The orthogonal projection from L'^{X, du) onto Tiy, is an integral operator K^r 
and H^i is a reproducing kernel Hilbert space of functions: 

$(x) = / K,{x,y)<l>{y)du{y), V<l> G H,. (24) 

The kernel is given explicitly by K„{x,y) =< i]„(^)\A^'^i]^(^y) >, if r]„(^y) G 
D{A'^^), My G X. So, the function $ G L^(X, dz/) is a wavelet transform 
(WT) iff it satisfies this reproducing relation. 

4. Reconstruction formula. 

The WT Wr^ may be inverted on its range by the adjoint operator, W~^ = W* 
on Ti^j to obtain for r]a{x) ^ D{A~^), \/x E X 

W-'<^ = f <l>(x)A;ir7.(,)dz/(x), $ G n,. (25) 

This is inverse WT. 

If A'^^ is bounded then So- is called a frame, if A^r = XI then S^ is called a tight 
frame. This two cases are generalization of a simple case, when 5*0- is an (ortho)basis. 

The most simple cases of this construction are: 
1. H = {e}. This is the standard construction of WT over a locally compact group. 
It should be noted that the square integrability of U is equivalent to U belonging 
to the discrete series. The most simple example is related to the affine {ax + b) 
group and yields the usual one-dimensional wavelet analysis 



TT 




{b,a)f]{x) = —f ^ . (26) 



For G = SIM{2) = R^ cxi (R+ x SO {2)), the similitude group of the plane, we 
have the corresponding two-dimensional wavelets. 

2. H = H^, the isotropy (up to a phase) subgroup of rj: this is the case of the 
Gilmore-Perelomov CS. Some cases of group G are: 
a). Semisimple groups, such as SU(N), SU(N|M), SU(p,q), Sp(N,R). 
b). the Weyl-Heisenberg group Gwh which leads to the Gabor functions, i.e. 
canonical (oscillator) coherent states associated with windowed Fourier transform 
or Gabor transform (see also part III A): 

K(g, P, v)f] (x) = exp{ifi{ip -p{x-q))f{x- q) (27) 

In this case H is the center of Gwh- In both cases time-frequency plane corresponds 
to the phase space of group representation. 

c). The similitude group SIM(n) of R"(n > 3): for H = SO{n — 1) we have the 
axisymmetric n-dimensional wavelets. 



d). Also we have the case of bigger group, containing both affine and Weyl-Hei- 
senberg group, which interpolate between affine wavelet analysis and windowed 
Fourier analysis: affine Weyl-Heisenberg group [11]. 

e). Relativity groups. In a nonrelativistic setup, the natural kinematical group 
is the (extended) Galilei group. Also we may adds independent space and time 
dilations and obtain affine Galilei group. If we restrict the dilations by the rela- 
tion oq = a^i where Oq, a are the time and space dilation we obtain the Galilei- 
Schrodinger group, invariance group of both Schrodinger and heat equations. We 
consider these examples in the next section. In the same way we may consider as 
kinematical group the Poincare group. When cq = a we have affine Poincare or 
Weyl-Poincare group. Some useful generalization of that affinization construction 
we consider for the case of hidden metaplectic structure in section III A. 

But the usual representation is not square-integrable and must be modified: 
restriction of the representation to a suitable quotient space of the group (the asso- 
ciated phase space in our case) restores square - integrability: G — > homogeneous 
space. 

Also, we have more general approach which allows to consider wavelets correspond- 
ing to more general groups and representations [12], [13]. 

Our goal is applications of these results to problems of Hamiltonian dynamics and 
as consequence we need to take into account symplectic nature of our dynamical 
problem. Also, the symplectic and wavelet structures must be consistent (this must 
be resemble the symplectic or Lie-Poisson integrator theory). We use the point of 
view of geometric quantization theory (orbit method) instead of harmonic analysis. 
Because of this we can consider (a) - (e) analogously. 

F Bases for Solutions 

We consider an important particular case of affine relativity group (relativity 
group combined with dilations) — affine Galilei group in n-dimensions. So, we have 
combination of Galilei group with independent space and time dilations: Gaff = 
Grn ixi -D2, where D2 = (R^)^ — R^, Gm is extended Galilei group corresponding 
to mass parameter m > {Gaff is noncentral extension of G txi D2 by R, where 
G is usual Galilei group). Generic element of Gaff is g = {^,bo,b]V]R,ao,a), 
where $ G R is the extension parameter in Gm, bo & H, b & R" are the time and 
space translations, v G R" is the boost parameter, R G SO{n) is a rotation and 
Co, a G R;^ are time and space dilations. The actions of g on space-time is then 
X H^ aRx + ttovt + b, t \-^ a^t + bo, where x = (xi, X2, ..., x„). The group law is 

a^ 1 

gg' = ($ H $' + avRb' + -aov'^b'^, bo + ao&o, b + aRb' + aovb'^; (28) 

Oo 2 

V H Rv' , RR'; aaa'o, aa) 

It should be noted that D2 acts nontrivially on Gm- Space-time wavelets associated 
to Gaff corresponds to unitary irreducible representation of spin zero. It may 



be obtained via orbit method. The Hilbert space is 7i = L^(R" x R, d/cdcj), 
k = {ki, ..., kn), where R" x R may be identified with usual Minkowski space and 
we have for representation: 

{U{g)'^){k,uj) = ^/a^expi{m^ + kb - ujbo)'^{k' ,uj'), (29) 

with k' = aR~^{k + mv), u' = aQ{uj — kv — ^mv'^), m' = {a^ /aojm. Mass m 
is a coordinate in the dual of the Lie algebra and these relations are a part of 
coadjoint action of Gaff- This representation is unitary and irreducible but not 
square integrable. So, we need to consider reduction to the corresponding quotients 
X = G/H. We consider the case in which H= {phase changes $ and space dilations 
a}. Then the space X = G/H is parametrized by points x = {bQ,b;v; R;aQ). 
There is a dense set of vectors r] E H admissible mod(//, ci/j), where ap is the 
corresponding section. We have a two-parameter family of functions /9(dilations): 
(3{x) = {no + X-faoY^"^, Ao,/io ^ R- Then any admissible vector r] generates a tight 
frame of Galilean wavelets 

Vmik,u;) = ^ao(/io + Aoao)"/V(^'-"'''^r7(A;',cu'), (30) 

with k' = {fXQ+ XoaY^"^ R^^ (k+mv) , uj' = aQ{uj—kv—mv'^/2). The simplest examples 
of admissible vectors (corresponding to usual Galilei case) are Gaussian vector: 
ri{k) ~ exp(— A;^/2m-u) and binomial vector: ri{k) ~ (1 + A;^/2m-u)~"/^, a > 1/2, 
where -u is a kind of internal energy. When we impose the relation oq = ^^ then 
we have the restriction to the Galilei-Schrodinger group Gs = Gm txi -Ds, where 
Ds is the one-dimensional subgroup of D2- Gg is a natural invariance group of 
both the Schrodinger equation and the heat equation. The restriction to Gg of the 
representation (29) splits into the direct sum of two irreducible ones U = U^®U^ 
corresponding to the decomposition L^(R" x R, dkduj) = 7i+ © 7Y_, where 

n± = L\D±,dkduj (31) 

= {^e L2(R" X R, dkdco), ^(A;, u) = ior u + k^ /2m = 0} 

These two subspaces are the analogues of usual Hardy spaces on R, i.e. the sub- 
spaces of ( ant i) progressive wavelets (see also below, part III A). The two repre- 
sentation U± are square integrable modulo the center. There is a dense set of 
admissible vectors f], and each of them generates a set of CS of Gilmore-Perelomov 
type. Typical wavelets of this kind are: 
the Schrodinger-Marr wavelet: 

7]{x, t) = {tdt + A)e"(-'+*')/2 (32) 



2m' 



the Schrodinger-Cauchy wavelet: 



So, in the same way we can construct invariant bases with explicit manifestation of 
underlying symmetry for solving Hamiltonian (12) or Lagrangian (19) equations. 



Ill SYMPLECTIC STRUCTURES, QUANTIZATION 
AND FAST WAVELET TRANSFORM 

A Metaplectic Group and Representations 

Let Spin) be symplectic group, Mp{n) be its unique two- fold covering - meta- 
plectic group [14]. Let V be a symplectic vector space with symplectic form ( , ), 
then R®V is nilpotent Lie algebra - Heisenberg algebra: 

[i?, V] = 0, [v, w] = (v, w) G R, [V, V] = R. 

Sp{V) is a group of automorphisms of Heisenberg algebra. 

Let N be a group with Lie algebra R(BV, i.e. Heisenberg group. By Stone- von 
Neumann theorem Heisenberg group has unique irreducible unitary representation 
in which 1 i— > i. Let us also consider the projective representation of simplectic 
group Spiy): Ug^Ug^ = c{gi,g2) ■ Ug-^g^, where c is a map: Sp{V) x SpiV) -^ S^, 
i.e. c is S'^-cocycle. 

But this representation is unitary representation of universal covering, i.e. meta- 
plectic group Mp{V). We give this representation without Stone- von Neumann 
theorem. Consider a new group F = A^' ixi Mp{V), cxi is semidirect product (we 
consider instead oi N = R®V the N' = S^ xV, S^ = {R/2txZ)). Let V* be dual 
to V, G{y*) be automorphism group of \/*.Then F is subgroup of G(y*), which 
consists of elements, which acts on V* by afSne transformations. 
This is the key point! 

Let gi, ..., g„;pi, ...,p„ be symplectic basis in V, a = pdq = J^PidQi and da be 
symplectic form on V*. Let M be fixed affine polarization, then for a & F the map 
a H^ ©a gives unitary representation of G: O^ : H{M) -^ H{M) 

Explicitly we have for representation of N on H(M): 

(6 J)*(x) = e-^-/(x), ejix) = fix - p) 

The representation of N on H(M) is irreducible. Let A^, Ap be infinitesimal opera- 
tors of this representation 

df 

then Agf{x)=i{qx)f{x), Apf{x) =^pj — {x) 

Now we give the representation of infinitesimal basic elements. Lie algebra of the 
group F is the algebra of all (nonhomogeneous) quadratic polynomials of (p,q) 
relatively Poisson bracket (PB). The basis of this algebra consists of elements 
l,qi,...,qn,Pu...,Pn,qiqj,qiPj,PiPj, ij = l,...,n, i<j, 



PB is {f:g} = ^-Q—Q---Q—Q- and {1,^} = for all g, 

{Pi, Qj} = ^ij, {PiQj, Qk} = SikQj, {Piqj,Pk} = -SjkPi, {PiPj,Pk} = 0, 
{PiPj, Qk} = SikPj + 5jkPi, {qiqj, Qk} = 0, {qiqj,Pk} = -^ikQj - SjkQi 

so, we have the representation of basic elements / i— * Aj : 1 t— > i, g^ i-^ ixk, 

This gives the structure of the Poisson manifolds to representation of any (nilpotent) 
algebra or in other words to continuous wavelet transform. 

The Segal-Bargman Representation. Let z = 1/a/2 ■ {p — iq), z = l/y^ ■ 
(p + iq), p = {pi, ■■■,Pn), Fn is the space of holomorphic functions of n complex 
variables with (/, /) < oo, where 



Consider a map U : H ^ Fn , where H is with real polarization, F„ is with complex 
polarization, then we have 

{U^){a)= I A{a,q)^'{q)dq, where A(a, g) = 7r-"/Vi/=^("'+^')+^"'? 

i.e. the Bargmann formula produce wavelets. We also have the representation of 
Heisenberg algebra on F„ : 

and also : uj = dj3 = dp A dq, where f3 = izdz. 

Orbital Theory for Wavelets. Let coadjoint action be < g ■ f,Y >=< 

f,Ad{g)^^Y >, where <, > is pairing g E G, f ^ g*, Y E Q. The orbit is 
Of = G ■ f = G/G{f). Also, let A=A(M) be algebra of functions, V(M) is A- 
niodule of vector fields, A^ is A-module of p-forms. Vector fields on orbit is 



a(O,X)X0) = ^(0(exptX/)) 



t=o 



where (p G A{0), f E O. Then Of are homogeneous symplectic manifolds with 
2-form n{(T{0, X)f, a{0, Y)f) =< f, [X, Y] >, and dQ = 0. PB on O have the next 
form {'^i,'^2} = p{'^i)'^2 where p is A^{0) -^ V{0) with definition Q{p{a),X) 
= i{X)a. Here ^i, ^2 e -4(C) and A{0) is Lie algebra with bracket {,}. Now let 
N be a Heisenberg group. Consider adjoint and coadjoint representations in some 
particular case. A^ = {z,t) E C x R, z = p+iq; compositions in N are {z, t) ■ (z', t') = 
{z + z',t + t' + B{z,z')), where B{z,z') =pq' — qp'. Inverse element is {—t,—z). Lie 



algebra n of N is (C, r) G C x i? with bracket [(C, r), {(', r')] = (0, B{C, C'))- Centre 
is 5 G n and generated by (0,1); Z is a subgroup exp z. Adjoint representation N on 
n is given by formula Ad{z, t){(, r) = (C, t + B{z, ()) Coadjoint: for / G n*, g = 
iz,t), (g- f)iCX) = /(C,^)- 5(^,0/(0,1) then orbits for which /jj ^ are plane 
in n* given by equation /(0, 1) = /i . If X = (C, 0), Y = {(', 0), X,Y e n then 
symplectic structure is 

n{a{0, X)f, a{0, Y)f) =< f, [X, Y] >= /(O, B{C, C))f^B{C, (') 

Also we have for orbit (9^ = N/Z and O^ is Hamiltonian G-space. 

According to this approach we can construct by using methods of geometric 
quantization theory many "symplectic wavelet constructions" with corresponding 
symplectic or Poisson structure on it. Very useful particular spline-wavelet ba- 
sis with uniform exponential control on stratified and nilpotent Lie groups was 
considered in [13]. 

B Applications to Melnikov Functions Approach 

We give now some point of applications of wavelet methods from the preced- 
ing parts to Melnikov approach in the theory of homoclinic chaos in perturbed 
Hamiltonian systems for examples from [1]. 

In Hamiltonian form we have: 

x = J-VH{x)+eg{x,e), & = oo, (x, 9) G i?" x T™, 

for e = we have: 

x = J-VH{x), = cj (34) 

For £ = we have homoclinic orbit XQ{t) to the hyperbolic fixed point Xq. For 
£ 7^ we have normally hyperbolic invariant torus T^ and condition on transver- 
sally intersection of stable and unstable manifolds W^{T^) and W^lT^) in terms of 
Melnikov functions M(0) for XQ{t): 

oo 

M(e) = / VH{xo{t)) A g{xo{t),ujt + Q)dt 

— oo 

This condition has the next form: 

M(eo) = o, ^,;^.J^M(eo)^o 

According to the approach of Birkhoff-Smale- Wiggins we determined the region in 
parameter space in which we can observe the chaotic behaviour [4]. 



If we cannot solve equations (34) explicitly in time, then we use the wavelet ap- 
proach from paper [1] for the computations of homoclinic (heteroclinic) loops as 
the wavelet solutions of system (34). For computations of quasiperiodic Melnikov 
functions 

f-niT 

M"/"(to) = / DH{x^{t)) A g{x^{t),t + to)rft 

JO 

we used periodization of wavelet construction from paper [1]. 

We also used symplectic Melnikov function approach in which we have: 



Mi{z) = \im. I {hi,h]^u,z)dt 



-T 



d,{z,e) = h,{z^) - K{zl) = eM,{z) + 0{e^) 

where {, } is the Poisson bracket, di{z,e) is the Melnikov distance. So, we need 
symplectic invariant wavelet expressions for Poisson brackets. The computations 
are produced according to invariant calculation of Poisson brackets, which is based 
on consideration in part III A and on operator representation from part III F (see 
below) . 



C Floer Approach for Closed Loops 

Now we consider the generalization of wavelet variational approach to the sym- 
plectic invariant calculation of closed loops in Hamiltonian systems [15]. As we 
demonstrated in [4] we have the parametrization of our solution by some reduced 
algebraical problem but in contrast to the cases from paper [1], where the solution 
is parametrized by construction based on scalar refinement equation, in symplec- 
tic case we have parametrization of the solution by matrix problems - Quadratic 
Mirror Filters equations. Now we consider a different approach. 

Let(M, u;) be a compact symplectic manifold of dimension 2n, u; is a closed 2- 
form (nondegenerate) on M which induces an isomorphism T*M ^ TM. Thus 
every smooth time-dependent Hamiltonian if : R x M ^ R corresponds to a 
time-dependent Hamiltonian vector field Xh : R x M — >■ TM defined by 

u{XHit,x),^) = -d,Hit,x)^ (35) 

for C^ G TxM. Let H (and Xh) is periodic in time: H(t + T,x) = H(t,x) and 
consider corresponding Hamiltonian differential equation on M: 

x{t) =XH{t,x{t)) (36) 

The solutions x{t) of (36) determine a 1-parameter family of diffeomorphisms ipt G 
Diff(M) satisfying iljt{x{0)) = x{t). These diffeomorphisms are symplectic: u = 



ip*uj. Let L = LrpM be the space of contractible loops in M which are represented 
by smooth curves 7 : R ^ M satisfying 7(t + T) = 7(t). Then the contractible T- 
periodic solutions of (36) can be characterized as the critical points of the functional 
= St '■ L — > IX'. 

Sril) = - f u*co+ r Hit, 7(t))dt, (37) 

Jd Jo 

where D C C be a closed unit disc and u : D ^ M is a smooth function, which on 
boundary agrees with 7, i.e. ■u(exp{27ri6}) = 7(6T). Because [u], the cohomology 
class of u, vanishes then S't(7) is independent of choice of u. Tangent space T^L is 
the space of vector fields ^ G C'^{'~f*TM) along 7 satisfying C,{t + T) = C,{t). Then 
we have for the 1-form d/ : TL -^ R 



dSrm = rHi, + dH{t, 7)0dt (3^ 

Jo 



and the critical points of S are contractible loops in L which satisfy the Hamiltonian 
equation (36). Thus the critical points are precisely the required T-periodic solution 
of (36). 

To describe the gradient of S we choose a on almost complex structure on M 
which is compatible with uj. This is an endomorphism J G C°°(End(TM)) satisfy- 
ing J^ = —I such that 

gi^,r])=uj{C,Jix)r]), C,V^T,.M (39) 

defines a Riemannian metric on M. The Hamiltonian vector field is then represented 
by Xnityx) = J{x)'VH(t,x), where V denotes the gradient w.r.t. the x-variable 
using the metric (39). Moreover the gradient of S w.r.t. the induced metric on L 
is given by 

grad5(7) = J(7)7 + Vif(t,7), 7 e L (40) 

Studying the critical points of S is confronted with the well-known difficulty that 
the variational integral is neither bounded from below nor from above. Moreover, 
at every possible critical point the Hessian of / has an infinite dimensional positive 
and an infinite dimensional negative subspaces, so the standard Morse theory is 
not applicable. The additional problem is that the gradient vector field on the loop 
space L 

d 

—7 = -grad/(7) (41) 

does not define a well posed Cauchy problem. But Floer [15] found a way to 
analyse the space A1 of bounded solutions consisting of the critical points together 
with their connecting orbits. He used a combination of variational approach and 



Gromov's elliptic technique. A gradient flow line of / is a smooth solution m : R ^ 
M of the partial differential equation 

- + J{u)- + VH{t,u)=0, (42) 

which satisfies u{s,t + T) = u{s,t). The key point is to consider (42) not as the 
flow on the loop space but as an elliptic boundary value problem. It should be 
noted that (42) is a generalization of equation for Gromov's pseudoholomorphic 
curves (correspond to the case VH = in (42)). Let Air = A4t{H, J) the space 
of bounded solutions of (42), i.e. the space of smooth functions u : C/iTZ -^ M, 
which are contractible, solve equation (42) and have finite energy flow: 

For every u G Mt there exists a pair x, y of contractible T-periodic solutions of 
(36), such that -u is a connecting orbit from y to x: 

lim u(s,t) =y(t), lim = x(t) (44) 

s— > — oo s— »+oo 

Then the approach from [1], which we may apply or on the level of standard bound- 
ary problem (42) or on the level of variational approach (43) and representation 
of operators (in our case, J and V) according to part III F (see below) lead us to 
wavelet representation of closed loops. 

D Quasiclassical Evolution 

Let us consider classical and quantum dynamics in phase space Q = R"^™- with 
coordinates {x,C,) and generated by Hamiltonian 7i{x,C,) G C°°(r2;i?). If ^J*^ : 
Q — > Q is (classical) flow then time evolution of any bounded classical observable 
or symbol b{x,^) e C°°{n,R) is given by bt{x,^) = 6(<I)J^(a;,0)- Let H = Op^iH) 
and B = Op^{b) are the self-adjoint operators or quantum observables in L'^{R"'), 
representing the Weyl quantization of the symbols 7i, b [14] 

where u G S{R^) and Bt = ^^Ht/nj^^-iHt/n -\^^ ^-^^q Heisenberg observable or quantum 
evolution of the observable B under unitary group generated by H . Bt solves the 
Heisenberg equation of motion 

i3t='-[H,Bt]. 
Let 6j(x, S,\ K) is a symbol of Bt then we have the following equation for it 



bt = {n,h}M, (45) 

with initial condition bo{x,C^, Ti) = b{x,C^). Here {f,g}M{x,i) is the Moyal brackets 
of the observables f,g e C°^(i?^"), {f,g}M{x,0 = fh - 5'tt/, where f]\g is the 
symbol of the operator product and is presented by the composition of the symbols 

f,g 

g{x + r,T + ^)dpdrdrda'. 

For our problems it is useful that {f,g}M admits the formal expansion in powers 
oih: 

{f,g}M{x,0 - {f,g} + 2'' E (-1)"" ■ id^fD^) ■ {d^gD'^J), 

|a+/3|=i>l 

where a = (ai, . . . , «„) is a multi-index, \a\ = ai + . . . + «„, D^ = —ihdx- So, 
evolution (45) for symbol bt{x,^]h) is 

bt = {n, b,} + ^ E (-1)'^' ■ ^'{d^nD%) ■ {dlhD'^H). (46) 

l"l+/3|=i>i 

At /i = this equation transforms to classical Liouville equation 

bt = {nA}- (47) 

Equation (46) plays a key role in many quantum (semiclassical) problem. We 
note only the problem of relation between quantum and classical evolutions or how 
long the evolution of the quantum observables is determined by the corresponding 
classical one [14]. Our approach to solution of systems (46), (47) is based on 
our technique from [l]-[7] and very useful linear parametrization for differential 
operators which we present in section III F. 

E SYMPLECTIC HILBERT SCALES VIA WAVELETS 

We can solve many important dynamical problems such that KAM perturba- 
tions, spread of energy to higher modes, weak turbulence, growths of solutions of 
Hamiltonian equations only if we consider scales of spaces instead of one functional 
space. For Hamiltonian system and their perturbations for which we need take into 
account underlying symplectic structure we need to consider symplectic scales of 
spaces. So, if u{t) = J'S/K{u{t)) is Hamiltonian equation we need wavelet descrip- 
tion of symplectic or quasicomplex structure on the level of functional spaces. It is 
very important that according to [16] Hilbert basis is in the same time a Darboux 
basis to corresponding symplectic structure. We need to provide Hilbert scale {Zg} 



with symplectic structure [16], [17]. All what we need is the following. J is a linear 
operator, J : Z^o -^ Z^o, J (Zoo) = Z^o, where Z^o = riZg. J determines an isomor- 
phism of scale {Zg} of order cij > 0. The operator J with domain of definition Z^o 
is antisymmetric in Z: < Jzi, Z2 >z= — < -^i, Jz2 >z, Z\, Zi G Z^o- Then the triple 

{Z, {Zs|s G -R}, a =< Jdz, dz >} 

is symplectic Hilbert scale. So, we may consider any dynamical Hamiltonian prob- 
lem on functional level. As an example, for KdV equation we have 

Z, = {uix) G H\T'^)\ I ^ u{x)<lx = 0}, s G R, J = d/dx, 

Jo 

J is isomorphism of the scale of order one, J = — ( J)~^ is isomorphism of order —1. 
According to [18] general functional spaces and scales of spaces such as Holder- 
Zygmund, Triebel-Lizorkin and Sobolev can be characterized through wavelet co- 
efficients or wavelet transforms. As a rule, the faster the wavelet coefficients decay, 
the more the analyzed function is regular [18]. Most important for us example is 
the scale of Sobolev spaces. Let Hk{R^) is the Hilbert space of all distributions 
with finite norm 

Let us consider wavelet transform 

WJ{b,a)= f dxlg f ^i^) f(x), 

b G -R", a > 0, w.r.t. analyzing wavelet g, which is strictly admissible, i.e. 

Cg,g = / — \g{ak)\^ < oo. 
Jo a 

Then there is a c > 1 such that 

dbda 



^llk(«") < /,„ ^(1 + a-'')\WAb,<^)\' < c\\s\\l^i 



This shows that localization of the wavelet coefficients at small scale is linked to 
local regularity. 

So, we need representation for differential operator (J in our case) in wavelet 
basis. We consider it in the next section. 

F FAST WAVELET TRANSFORM FOR 
DIFFERENTIAL OPERATORS 

Let us consider multiresolution representation . . . C V2 C Vi C Vq C K_i C 
V-2 ■ ■ ■ (see our other paper from this proceedings for details of wavelet machinery). 



Let T be an operator T : L^{R) —^ L^{R), with the kernel K{x, y) and Pj : L'^{R) — * 
Vj {j G Z) is projection operators on the subspace Vj corresponding to j level of 
resolution: 

k 

Let Qj = Pj-i — Pj is the projection operator on the subspace Wj then we have the 
following "microscopic or telescopic" representation of operator T which takes into 
account contributions from each level of resolution from different scales starting 
with coarsest and ending to finest scales: 

T = J2(QjTQj + QjTP, + PjTQj). 

We remember that this is a result of presence of affine group inside this construction. 
The non-standard form of operator representation [19] is a representation of an 
operator T as a chain of triples T = {Aj,Bj,rj}jfzz, acting on the subspaces Vj 
and Wj-. 

A, : W, -. W,,B, : V, ^ W„T, : W, ^ \/„ 

where operators {Aj, Bj, Tj}j(zz are defined as 

The operator T admits a recursive definition via 

T = ( "^^+1 ^^+1 I 

where Tj = PjTPj and Tj works on Vj : Vj ^ Vj. It should be noted that operator 
Aj describes interaction on the scale j independently from other scales, operators 
BjjTj describe interaction between the scale j and all coarser scales, the operator 
Tj is an "averaged" version of Tj_i. 

The operators Aj, Bj, Tj, Tj are represented by matrices a\ [3^ , 7-', s^ 



(^i,k' = J J K{x,y)-^j,k{x)-4jj,k'{y)dxdy 

Pk,k' = J J K{x, y)^j,k{x)(Pj,k'{y)dxdy (4J 

lk,k' = J J K{x,y)ipj,kix)^ljj,k'iy)dxdY 

4,k'= / K{x,y)^j^k{x)^j,k'{y)dxdy 



We may compute the non-standard representations of operator d/dx in the wavelet 
bases by solving a small system of linear algebraical equations. So, we have for 
objects (48) 



(3l^ = 2-^ [■^{2-^x - i)ip'{2-^x - £)2-Mx = 2^^-^ 
^i^ = 2-^ f ip{2-^x - i)i)'{2-^x - £)2'Mx = 2-J7i^^, 



where 



Oi£= ^p(x — i)-—^p(x)dx 
dx 



(3e= I iIj{x - i)—ip{x)dx 



7£ = / ip{x - i)—ilj{x)dx 
then by using refinement equations 

L-l 

ip{x) = V^Y^ hk(p{2x - k), 

L-l 

^{x) = ^/2j2 fl'fcV^(2x - k), 



k=0 



9k 



-D'^h 



L-k-\i 



where 



A; = 0, . . . , L — 1 we have in terms of fihers (/i^, Qk)'- 

L-l L-l 

ttj = 2 ^ ^ gkgk'r2i+k-k', 

k=0 k'=0 
L-l L-l 

/5j = 2 2^ 2^ gkhk'r2i+k-k', 

k=0 k'=0 
L-l L-l 

7i = 2 XI 51 hkgk'r2i+k-k', 

k=0 k'=0 



re = I ip{x — £) —ip{x)dx, £ G Z. 



Therefore, the representation of d/dx is completely determined by the coefficients 
r£ or by representation of d/dx only on the subspace Vq. The coefficients ri,i & Z 
satisfy the following system of linear algebraical equations 



re 



1 



L/2 



^2« + - XI (^2k-i{r2i-2k+i + r2e+2k-i) 
^ k=l 



and J2e^re = —1, where a2k-i = 2 J2i=o'' hihi+2k-i, k = 1, . . . , L/2 are the autocor- 
relation coefficients of the filter H. If a number of vanishing moments M > 2 then 



this linear system of equations has a unique solution with finite number of non-zero 
r£, r£ ^ for — L + 2 < i < L — 2,r£ = —r^i. For the representation of operator 
dP" /dx^ we have the similar reduced linear system of equations. Then finally we 
have for action of operator Tj{Tj : Vj — > Vj) on sufficiently smooth function /: 

k&z \ e / 

where ^j,k{x) = 2~^/'^(p(2~^x — k) is wavelet basis and 



are wavelet coefficients. So, we have simple linear parametrization of matrix rep- 
resentation of our differential operator in wavelet basis and of the action of this 
operator on arbitrary vector in our functional space. Then we may use such repre- 
sentation in all preceding sections. 

IV MAPS AND WAVELET STRUCTURES 

A Veselov-Marsden Discretization 

Discrete variational principles lead to evolution dynamics analogous to the Euler- 
Lagrange equations [9]. Let Q be a configuration space, then a discrete Lagrangian 
is a map L : Q x Q — * R. usually L is obtained by approximating the given 
Lagrangian. For N G N^ the action sum is the map S : Q^^^ — > R defined by 

N-l 

S= J2L{qk+i,qk), (49) 

fc=0 

where qk & Q, k > 0. The action sum is the discrete analog of the action integral 
in continuous case. Extremizing S over gi, ...,gAr_i with fixing qo,qN we have the 
discrete Euler-Lagrange equations (DEL): 



^2^(gfc+i, Qk) + DiiQk, Qq-i) = 0, (50) 



for A; = l,...,iV-l. 
Let 



u, ..., . 



<l>:QxQ^QxQ (51) 

and 

^{qk,qk~i) = {qk+i,qk) (52) 

is a discrete function (map), then we have for DEL: 



D2L o $ + DiL = (53) 

or in coordinates g* on Q we have DEL 

^— o$(gfc+i,gfc) + — - — {qk+i,qk) = 0. (54) 

It is very important that the map $ exactly preserves the symplectic form u: 

^ = ainj (^fc+1, gfe)dgfc A dg^+1 (55) 



B Generalized Wavelet Approach 

Our approach to solutions of equations (54) is based on applications of general 
and very efficient methods developed by A. Harten [20], who produced a "General 
Framework" for multiresolution representation of discrete data. It is based on con- 
sideration of basic operators, decimation and prediction, which connect adjacent 
resolution levels. These operators are constructed from two basic blocks: the dis- 
cretization and reconstruction operators. The former obtains discrete information 
from a given continuous functions (flows), and the latter produces an approxima- 
tion to those functions, from discrete values, in the same function space to which 
the original function belongs. A "new scale" is defined as the information on a 
given resolution level which cannot be predicted from discrete information at lower 
levels. If the discretization and reconstruction are local operators, the concept of 
"new scale" is also local. The scale coefficients are directly related to the prediction 
errors, and thus to the reconstruction procedure. If scale coefficients are small at a 
certain location on a given scale, it means that the reconstruction procedure on that 
scale gives a proper approximation of the original function at that particular loca- 
tion. This approach may be considered as some generalization of standard wavelet 
analysis approach. It allows to consider multiresolution decomposition when usual 
approach is impossible ((5-functions case). 

Let F be a linear space of mappings 

Fc{/|/:X^r}, (56) 

where X, Y are linear spaces. Let also D^ be a linear operator 

A:/^^'}, v^ = Dkf, v'^ = {v1}, v^eY. (57) 

This sequence corresponds to k level discretization of X. Let 

Du{F) = V^ = spanJT^f } (58) 

and the coordinates of v^ G V^ in this basis are v^ = {'Of}, v^ & S'': 



V 



k \ ' X'.k^k 



E^'vt (59) 



Dk is a discretization operator. Main goal is to design a multiresolution scheme 
(MR) [20] that apphes to all sequences s E S^, but corresponds for those sequences 
^L ^ ^L^ which are obtained by the discretization (56). 

Since Dk maps F onto V'' then for any v'' C V'' there is at least one f in F such 
that Dkf = v^. Such correspondence from / G -F to f '^ G l^'^ is reconstruction and 
the corresponding operator is the reconstruction operator R^: 

Rk-.Vk^ F, DkRk = Ik, (60) 

where Ik is the identity operator in V^ {R^ is right inverse of D^ in V'^). 

Given a sequence of discretization {Dk} and sequence of the corresponding re- 

>k 
k-l 



construction operators {Rk}-, we define the operators D^ ^ and Pf" 



D^l = Dk-iRk ■■ Vk -^ Vk-i (61) 

-Pfe-l = DkRk-l '■ Vk-1 -^ Vk 

If the set Dk in nested [20] , then 

Dl-'Pt^ = 4-1 (62) 

and we have for any f & F and any p G -F for which the reconstruction Rk-i is 
exact: 

Dl-\Dkf) = Dk-J (63) 

Pt,{Dk-ip) = DkP 

Let us consider any v^ G V^, Then there is / G -F such that 

v"- = DlI, (64) 

and it follows from (63) that the process of successive decimation [20] 

^k-i ^jjk-i^k^ A; = L,...,1 (65) 

yields for all k 

v^ = Dkf (66) 

Thus the problem of prediction, which is associated with the corresponding MR 
scheme, can be stated as a problem of approximation: knowing Dk-if, f E F, find 
a "good approximation" for Dkf. It is very important that each space V^ has a 
multiresolution basis 

^M = {0°'n.{{#%}ti (67) 



and that any v'^ G V^ can be written as 

^;^ = E«°0°•VEE4#^ (68) 

i fc=l j 

where {d''} are the k scale coefficients of the associated MR, {v^} is defined by (59) 
with A; = 0. If {-Dfc} is a nested sequence of discretization [20] and {Rk} is any 
corresponding sequence of hnear reconstruction operators, then we have from (68) 
for v^ = Dif applying Rl. 

RlDj = E M^' + E E 4<'^ (69) 

i fc=l j 

where 

0- '"^ = i?L0°'^ e F, ^)'^ = i?^^^^'^ G F, /^o/ = E/^°- (70) 

When L ^ oo we have sufficient conditions which ensure that the limiting process 
L — i> cxD in (69, 70) yields a multiresolution basis for F. Then, according to (67), 
(68) we have very useful representation for solutions of equations (54) or for different 
maps construction in the form which are a counterparts for discrete (difference) 
cases of constructions from paper [1] . 
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